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Abstract 


The Hybrid Wing Body (HWB) aircraft is of great interest for future transport concepts due to its 
promises of reduced aircraft noise, nitrous-oxide emissions, and fuel consumption. A design parameterization 
method for HWB configurations with mail slot nacelle has been developed for a fast exploration of design 
space in conceptual and preliminary design phases of a HWB configuration. A HWB planform model by 
Laughlin" was implemented, and the Class Shape Transformation (CST) airfoil generation method by 
Kulfan”” was utilized to construct the needed geometry for computational high fidelity aerodynamic 
simulations. Geometric constraints for the parameterization such as internal cabin and cargo hold layouts 
were imposed on the geometry generation. A CFD simulation was performed for a HWB configuration 
generated by the current geometric modeler, clearly showing a significant effect of the installed nacelle on the 
flowfield. 
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Nomenclature 


Nondimensional CST z coordinates for airfoil upper surface 

Nondimensional CST z coordinates for airfoil lower surface 

Nondimensional y coordinates 

Nondimensional local y coordinate 

Sweep angle of the leading edge of the planform 

Sweep angle of the outboard wing of the planform 

Coefficient representing the distance between the cabin body and the root chord of the outboard wing 
Coefficient representing the reach of the curved section of the leading edge of the outboard wing 
Coefficient representing the reach of the curved section of the trailing edge of the outboard wing 
Nondimensional trailing edge thickness for the upper surface of the nth control airfoil on the aircraft 
Nondimensional trailing edge thickness for the lower surface of the nth control airfoil on the aircraft 
Nondimensional CST x coordinates 

CST optimization variable vector for the upper surface of the airfoil 

CST optimization variable vector for the lower surface of the airfoil 
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bow = Span of the outboard wing 


b = Total span 

Clocal = Local chord length 

Crcb = Root chord of the cabin body 

Crew = Root chord of the outboard wing 

Cheb = Tip chord of the cabin body 

Chow = Tip chord of the outboard wing 

h = Offset vector of the control airfoils for the aircraft 

Anacelle = Height of the nacelle as referenced from the surface of the aircraft 
lew = Length of the cabin body 

N = Bernstein polynomial order for CST airfoil generation of the aircraft 
Ni = Class function exponent for CST airfoil generation of the aircraft 
N> = Class function exponent for CST airfoil generation of the aircraft 
Web = Width of the cabin body 

XLE = Leading edge position vector 

xX = Dimensionalized CST x coordinates 

Y = Dimensionalized CST y coordinates 

Lup = Dimensionalized CST z coordinates for airfoil upper surface 

Ligyi = Dimensionalized CST z coordinates for airfoil lower surface 


I. Introduction 


Hybrid Wing Body (HWB) aircraft is characterized by a flattened and airfoil-shaped body, which produces a 
substantial portion of the total lift. The body form is composed of distinct and separate wing structures, though the 
wings are smoothly blended into the body. This concept has been studied widely and results suggest remarkable 
performance improvements over the conventional tube and wing transport!”. HWB incorporates design features 
from both a futuristic fuselage and flying wing design, which houses most of the crew, payload and equipment 
inside the main centerbody structure. The purported advantages of the HWB approach are efficient high-lift wings 
and a wide-shaped body, which together reduce the wing loading and boost the spanwise lift distribution. This 
enables the entire craft to contribute to lift generation with the result of potentially increased fuel efficiency and 
range; the smoothly blended wing-body reduces interference drag. In addition, the wide airfoil-shaped body has 
demonstrated shielding downward propagating noise from the engines mounted on the upper surface of the vehicle’. 
Lyu et al. pointed out the challenges in designing HWB and reviewed major studies done in design optimization of 
the HWB configuration. Those efforts are very important in designing the unconventional HWB due to the fact that 
little historical data can be drawn, unlike the conventional tube-and-wing transport aircraft. However, the embedded 
engines in HWB vehicle, instead of pylon-mounted engine, would be partly submerged in boundary layer from 
center body upstream of the inlet with boundary-layer-ingestion (BLI) offset inlet employed. The BLI offset inlet 
provides better propulsion efficiency by reducing form drag, wetted area, structural weight and noise. However, the 
higher flow distortion and lower pressure recovery at the Aerodynamic Interface Plane (AIP) found in the BLI offset 
inlet are disadvantages. Consequently the low momentum flow at the low recovery region within the highly 
distorted flows may separate inside the diffuser and results in insufficient air supply to fan face®. Moreover, the 
installation of the nacelle is shown to significantly degrade the potential of the highly desirable aerodynamic 
efficiency of HWB in clean wing, see Ref. [7]. Hence, the main premise of the present work is to develop a design 
approach that includes both the airframe and propulsion systems at the outset of aerodynamic design of the aircraft. 
A similar propulsion-airframe integration (PAI) design approach was conducted on an over-wing nacelle airliner 
configuration®. This paper will illustrate an efficient geometric parameterization of a complex propulsion-airframe 
integration of the HWB, using the N3-X*"" hybrid wing-body with a mail-slot inlet/nozzle nacelle as a test case. 


To determine the optimal aerodynamic performance of the aircraft, a geometric model represented by a set of 
parameters (design variables) must be created in which these design variables can be quickly and easily altered. The 
Class Shape Transformation (CST) '', a non-dimensional airfoil generation method, is utilized to construct the 
needed geometry for preliminary aerodynamic testing. The remainder of this paper is organized as follows. Section 
II illustrates the geometry definition of HWB planform and its associated parameters by using CST. Section II 
presents our geometric representation of HWB clean wing body and nacelle, including the validation with the N3-X 
clean wing body. Section IV gives the demonstration of aerodynamic analysis from the computational fluid 


dynamics (CFD) simulations. In this study, we employ GO-flow’, a high-fidelity Euler / Navier-Stokes flow solver. 
Finally, concluding remarks are given in Section V. 


II. HWB Planform Overview 


The planform is the two-dimensional shadow of the aircraft as viewed from a position directly above the craft. It 
is an ideal starting point for the design because it serves as the two-dimensional basis by which the CST method can 
be employed to generate the three-dimensional aircraft body. The following HWB planform design was partially 
derived from the work of Laughlin’s'' with the HWB planform. 


A. Definition of Planform Shape 


Shown in Fig. 1 is the HWB planform with all control points labeled. The inboard sections of the aircraft 
planform are defined via splines passing through the inboard control points and the outboard sections are defined 
linearly. Here, a cubic spline for [y1, y3], through the first three (control) points and a linear spline through the last 
two are created. Details for the determination of the leading edge and trailing edge splines are given in Appendix A. 
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Figure 1: Planform definition of Hybrid Wing Body configuration 


B. Planform Control Points 


The HWB aircraft planform is composed of ten control points whose locations are derived from the twelve 
planform design variables: Cycb, Cicbs Crows Cows Webs lebs Dow, Ai, Az, Ai, Az, and A3. These ten control points, (x),y1) to 
(X10,¥10), are found through solving the following two systems of linear equations: Eqs. (1)-(10) for the x coordinates 
and Eqs. (11)-(20) for the y coordinates. 
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WI. The Clean Wing HWB Aircraft and Nacelle 
III-1 : The Clean Wing HWB Aircraft 


The first step in the generation of the three-dimensional HWB aircraft is the creation of the clean airframe. This 
is the body of the aircraft sans any additional external features, such as nacelles. To generate the clean wing body, 
four control sectional airfoils are placed along the spanwise direction at yl, y3, y7 and y9, respectively. 

The airfoil shape at each design section is generated by the CST method. Although the CST method can be used 
either in one-dimensional formulation for airfoils or two-dimensional formulation for wing shapes, in this study the 
design section airfoils are defined by the two-dimensional CST formulation, and the sectional shapes between the 
design sections are generated by interpolation of CST parameters at the two design sections. This is because the 
HWB configuration has abrupt change in sectional shapes in the spanwise direction that is hard to capture with only 
four design sections. 


A. Control Airfoils 


The control airfoils are the elements within the HWB aircraft design that determine the three-dimensional shape 
of the aircraft. They are created through the employment of the CST method and dimensionalized based on planform 
locational parameters. When defining the nondimensional airfoil shape, a nondimensional local x coordinate, w, is 
created; w is identical for each control airfoil and is defined as a nondimensional chord line that is equal to zero at 
the leading edge of the airfoil and one at the trailing edge. Once yw is defined, the CST method can be employed to 
determine the nondimensional z coordinates, Cup and Gow. Gup and Cow are defined for the nth control airfoil by Eqs. 
(21) and (22), representing the CST method. 
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In Eqs. (21) and (22), A,,, and A), refer to the i" value of the design parameter vector for the upper and lower 
surfaces of the section airfoil, and €,, and €, refer to the nondimensionalized twist. Here &,, = €, so that sectional 
airfoils have zero trailing edge thickness. 

To dimensionalize each sectional airfoil shape, the local chord length is used for x and z coordinates and the total 
span is used for the y coordinates. The y term is not used in airfoil generation; it is used solely to determine the 
spanwise location of each airfoil. Furthermore, the x and z coordinates each have an offset term added to them. The 
x coordinates have an offset term corresponding to the leading edge position of each respective CAF; this offset is a 
function of the planform. The z coordinates have an offset that is used to counteract the fact that the CST method 
defines the leading edge z coordinate at z = 0. This z offset can be utilized to create such designs as a dihedral or 
anhedral wing aircraft. The dimensionalizing procedure is shown in Eqs. (23)-(26). n is defined as 0 at the 
symmetric section and | at the wing tip. 


X= WCtocal + XLE (23) 
Y = 0.5nb (24) 
Lup ae Sup Clocal +h (25) 


Ziow = Stow€tocal +h (26) 


In this study, the number of design parameters is set as seven, thus N = 6 in Eqns. (21) and (22). As the initial 
sectional airfoil shape in this study, the sectional shapes of the N3-X configuration [9] is extracted and fitted by the 
CST method by gradient-based optimization. Figure 2 shows the results for the curve fitting for airfoils at sections y 
=0, 0.3, 0.7, and 1.0 of the N3-X airframe. 
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Figure 2: CST Control Airfoils for the Hybrid Wing Body Aircraft 


B. Interpolation of sectional airfoil shapes 


At any given spanwise location, the sectional airfoil must have its own set of unique CST parameters, twist, and 
z-offset values. A linear interpolation of the design parameters between two design sections was adopted for 
efficient determination of the parameters. However, it is important to note that the interpolated sectional airfoils are 
not simply linear interpolations of the design sectional airfoils. Since each airfoil is dimensionalized based on 
planform locational variables such as local chord length and spanwise position, the interpolated sectional airfoils 
create a smooth, planform-dependent, blended transition between any two design sectional airfoils. 

The CST parameters are determined via interpolation based on a local non-dimensional y coordinate, Nioca. When 
doing this interpolation in the middle of two design sections, the design section closer to the root chord is defined as 
Niocal = 0 and the design section closer to the tip is defined as Nioca) = 1. This definition provides for an extremely 
simple linear interpolation of design variables, as shown in Eqs. (27)-(31). 


€. = Sout (Sia Sou) Misca (27) 
& = Sor + (S10 — S01) Mocat (28) 
h = ho + (hy — ho) Mocat (29) 
Ay = Agu + (Aru - Aou)Mocal (30) 
Ay = Aga + (411 — 40,1) Mocat (31) 


Once the design parameters are obtained they are passed to Eqs. (21) and (22) to find the nondimensional z 
coordinates for each section airfoil. Following this, Eqs. (23)-(26) are utilized to dimensionalize each section airfoil. 
A front view and an isoperimetric view of the clean wing HWB aircraft is shown in Fig. 3 and Fig. 4, respectively. 


Figure 3: Front View of the Clean Wing Hybrid Wingbody Aircraft 


Figure 4: Isoperimetric View of the Clean Wing Hybrid Wing-body Aircraft 


IlI-2 Mail Slot Nacelle 


The final step in generating the geometry of the HWB aircraft is to create the nacelle. The nacelle for the N3-X 
HWB aircraft is dubbed the “mail-slot” nacelle. It is located on the upper surface of the aircraft (thus aiding in noise 
reduction because the actual body of the aircraft creates a shield for the noise) and houses the propulsion system. It 
was given the name “mail-slot” nacelle due to the fact that actual nacelle structure bears resemblance to a mail-slot. 


A. Nacelle Geometry 


The nacelle geometry was first generated with a straight leading edge and trailing edge. Translational 
movements can be made in the chordwise direction for any given nacelle planform on the airframe. Also, rotational 
movements for each sectional airfoil of the nacelle is required so the nacelle is aligned with the local slope of the 
airframe upper surface at the installation location. In this study, it was assumed that the mail slot nacelle has 
constant chord length along the span of the nacelle, and the nacelle trailing edge follows the trailing edge of the 
airframe. This assumption can be easily changed when needed. 

This procedure is the same as our previous work [9] for mail slot nacelle geometry generation with all the flow 
passages inside the nacelle. The difference is that the current nacelle has a simplified geometry with a single flow 
passage for the mail slot instead of having many number of internal flow passages or center body for propulsive fan 
installation. This simplified nacelle geometry is for fast turnaround time for geometry modeling and high fidelity 
flow simulations at conceptual and preliminary design phases of the HWB configuration with mail-slot nacelle, 
while keeping aerodynamic effects of the existence of mail slot on the airframe. The single passage is not a flow- 
through passage, but terminated by a boundary plane for numerical boundary conditions of specifying back pressure 
for a subsonic exit flow condition. The exit of the nacelle has another single plane for nozzle exit boundary 
condition as shown in Fig. 5. 


(b) Simplified mail-slot nacelle Structure 


Figure 5: Mail-slot nacelle geometry 


B. Formulating the Mail-slot Shape 


The three dimensional shape of the nacelle is mail-slot shaped, as noted above; the nacelle planform sits at a 
height hyacetie above the surface of the HWB aircraft. Therefore, the height of the nacelle is a function of the upper 
surface of the aircraft. Due to this, legs and rounded edges connecting the legs to the planform must be defined. The 
rounded edges are created as quarter ellipses in the y-z plane. Thus the nacelle legs have a height that is equivalent 
to the difference between the planform height at the nacelle edges and the height of the ellipse. Fig. 5 shows the 
mail-slot shape of the nacelle. 

Figure 6 shows a front view of the mail slot nacelle installed on the airframe. The figure clearly shows the 
horizontal region, quarter ellipse region and vertical leg region of the nacelle. 


Figure 6: Front view of mail slot nacelle: horizontal, quarter elliptic, and vertical leg regions of the nacelle 
cowl are clearly seen. 


C. Nacelle Airfoil Sections 


To promote aerodynamic efficiency, the cross section of the mail slot nacelle is also airfoil shaped. Therefore, 
the cross section of the nacelle can be created exactly as the control airfoils were created. The primary difference is 
that the nacelle airfoils are generated by three-dimensional formulation [10] of CST because of smooth variation of 
sectional shapes of the mail slot nacelle. 


IV. Geometric Constraints for HWB Airframe 


The airframe outer mold line (OML) defined by the present HWB planform and sectional shape parameterization 
has to be defined so that the airframe wraps around the internal cabin and cargo hold without any intersection with 
the internal volumetric structures. In order to specify the geometric constraints on the airframe configuration, the 
HWB cabin layout suggested by Nickol [12] for 301 passengers shown in Fig. 7 was adopted here. Based on the 
layout of the cabin design and height of the cabin suggested by Nickol, we generated a volumetric model of cabin 
and cargo hold as shown in Fig. 8. Design parameters such as the airframe planform, sectional airfoil thicknesses, 
and even the dimensions of the cargo hold were adjusted so that there is no intersection between the OML and the 
internal frames. The volume of the cargo hold was kept larger than the required volume of 5656 ft*, which is the 
cargo hold volume of Boeing 777-200LR. The root chord length of the HWB configuration is 118 ft, and the leading 
edge sweep angle of centerbody is set as 51 degrees following the cabin sweep angle in Fig. 7(a). As mentioned 
earlier, the sectional shapes are the same as the N3-X configuration [9] as shown in Fig. 2, and the maximum 
thicknesses of the design sections are used as additional design parameter. 

The HWB configuration in Fig. 8 is not any result of shape design optimization, but just a demonstration of the 
geometric parameterization method of the present study. 


(a) Top view of cabin layout 
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(b) Front view of internal layout of cabin and cargo hold 


Figure 7: HWB Cabin layout design for 301 passengers [12] 
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Figure 8: HWB configuration and cabin/cargo hold positioned in the airframe 


V. Demonstration of CFD Analysis 


The geometry model generated by the present HWB parameterization method can be used as a baseline 
geometry for unstructured surface and volume meshing for CFD simulations. Computational surface meshes with 
and without a nacelle is shown Figure 9. Once a computational mesh is generated, it can be deformed for another set 
of geometric parameters for a fast turnaround time of flow simulation during design optimization process of the 
HWB airplane. 

Go-flow, an inviscid/viscous three dimensional unstructured mesh flow solver, was adopted for flow 
simulations. In this study, inviscid simulations for Mach 0.84 and angle-of-attack of 3 degrees were conducted for a 
HWB configuration with and without nacelle on it. 
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Figure 9: Computationa meshe 
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(a) HWB configuration with nacelle 


(b) HWB configuration without nacelle 


Figure 10: Computed pressure contours for HWB with and without nacelle 


Figure 10 shows surface pressure contours and shock wave visualizations for the two configurations. There are 
shock waves on the outboard upper surface of the airframes. The outboard shock wave extends further inboard for 
the HWB without nacelle. The HWB with the nacelle has shock wave on the outer cowl surface which was also seen 
in our previous studies for the N3-X configuration. [9,13] 


VI. Conclusions 


A design parameterization method for HWB configurations with mail slot nacelle has been developed for a fast 
exploration of design space in conceptual and preliminary design phases of a HWB configuration. Geometric 


constraints for the parameterization such as internal cabin and cargo hold layouts are imposed on the geometry 
generation. A CFD simulation was performed to demonstrate the analysis for a HWB configuration generated by the 
current geometric modeler. The results clearly show the significant effect of the nacelle on the flowfield. For future 
work, the geometric modeling method will be applied to the conceptual design of a HWB configuration with 
turboelectric distributed propulsion system. 
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Appendix A 


We applied a piecewise splining algorithm specifically suited for this aircraft. The algorithm splines through given 
four control points T;(x;,y;), i= 1, 2, 3, 4 in the x-y plane. 


Figure A1: A piecewise splining algorithm through four points 


In the derivation, y will be presented as a function of x. In the planform and airfoil loft of the above, x was a 
function of y. Simply input x from the aircraft into y for the algorithm, and vice versa, to solve for the splines. 


The function f (x) is composed of three polynomials, of the format described below. 


SX) =a +), 0 =x) Poe x4)? d(x = ai)? (Al) 
Sp (x) = Ag + by (x — Xz) + C2 (x — im)? +a30= x2)8 (A2) 
S3(¢) = mx +k (A3) 


Spline S; connects points T; and T;,,. 


The third spline can be readily solved from T3 and T,. 


Pe Ya-Y3 (A4) 
X4-X3 
k =y, —Mx4 (A5) 


The first two splines are cubic piecewise interpolating splines with clamped end conditions. There are eight 
unknowns in Eqs. (A1) and (A2), so eight rules will be used to define them. It is convenient to use hy = Xj41 — Xj. 


1. The slope of spline one must be zero at x, 


S1'(%1) = b, = 0 (A6) 
2. The slope of spline two at x3 must equal the slope of spline three. 

S53 (x3) = by + 2c,h, + 3d,h2 =m = S3(x3) (A7) 
3. The slope of spline two must equal the slope of spline one at xz 

Sy (Xz) = Si (x2) (A8) 

by = by + 2c,h, + 3d,h? (A9) 

b, + 2c,h, + 3d,h? —b, =0 (A10) 
4. The second derivate of spline two must equal the second derivative of spline one at x, 

53’ (X2) = Sy! (x2) (Al]) 

2Cy = 2c, + 6d,h, (Al2) 

2c, + 6d,h, — 2c, = 0 (A13) 
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. The value of spline one at x, must equal y, 


51%) =a, =, (Al4) 
6. The value of spline one at xz must equal y, 


Si (%_) = a, t+ byh, +c,h? + dh? = y, (A15) 


7. The value of spline two at x. must equal yz 


52 (x2) = az = y2 (Al6) 
. The value of spline two at x3 must equal y3 
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So(x3) = a, + byh, +.c,h3 + d,h3 =y, (A17) 


These equations are all correctly formatted to create a linear system of equations. 
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Solving this system of equations will yield the correct coefficients for the cubic piecewise interpolating splines. The 
complete spline function follows rules that conform well to the outer curvature of the aircraft. 


It was found that the above algorithm results in non-realistic curves for trailing edges when the spanwise coordinates 
of control points T2 and T3 are close each other. For a more robust shape design, the order of the first polynomial 
S1 for the trailing edge spline curve was reduced from cubic to quadratic by setting the condition 4 as 


d, =0, (A19) 


and modifying the matrix equation (A18) accordingly. Figure A2 compares the effect of modified spline method on 
the resulting trailing edge curve. 
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Figure A2: Unrealistic trailing edge curve (left) obtained when control point P3 is located too close to P2 in 
the spanwise coordinate (left figure: cubic function for S,, right figure: quadratic function for S,) 


